Laser cooling to quantum degeneracy 
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We report on Bose-Einstein condensation (BEC) in a gas of strontium atoms, using laser cooling 
as the only cooling mechanism. The condensate is formed within a sample that is continuously 
Doppler cooled to below 1 /iK on a narrow-linewidth transition. The critical phase-space density 
for BEC is reached in a central region of the sample, in which atoms are rendered transparent for 
laser cooling photons. The density in this region is enhanced by an additional dipole trap potential. 
Thermal equilibrium between the gas in this central region and the surrounding laser cooled part 
of the cloud is established by elastic collisions. Condensates of up to 10 atoms can be repeatedly 
formed on a timescale of 100 ms, with prospects for the generation of a continuous atom laser. 



PACS numbers: 03.75. Jk, 37.10.De 

Laser cooling has revolutionized contemporary atomic 
and molecular physics in many respects, for example 
pushing the precision of clocks by orders of magnitude, 
and enabling ion quantum computation [1 . Since the 
early days of laser cooling, the question has been asked if 
the quantum degenerate regime could be reached using 
this efficient method as the only cooling process. De- 
spite significant experimental and theoretical effort to 
overcome the limitations of laser cooling this goal has 
been elusive. So far, laser cooling had to be followed by 
evaporative cooling to reach quantum degeneracy [2] . 

A gas of bosonic atoms with number density n and 
temperature T enters the quantum-degenerate regime 
and forms a Bose-Einstein condensate if its phase-space 
density n\^ B exceeds a critical value of 2.612. Here, 
Kib = /i/(27rmfcsT 1 ) 1 / 2 is the thermal de Broglie wave- 
length, where h and ks are Planck's and Boltzmann's 
constant, respectively, and m is the mass of an atom. 
Since n\^ B oc rtT~ 3 / 2 , low temperatures in combination 
with high densities have to be reached to obtain quan- 
tum degeneracy. Numerous studies, mainly carried out 
in the 1980's and 90's, have paved the ground to the 
present state of the art of laser cooling and have identi- 
fied the limitations of this technique [5] . 

The longstanding goal of reaching the quantum de- 
generate regime by laser cooling [IHS] can be discussed 
in terms of three main experimental challenges. First, 
temperatures in the low microkelvin regime have to be 
reached. Only here, quantum degeneracy can be ob- 
tained at a density that is low enough to avoid fast de- 
cay of the gas by molecule formation. This challenge has 
been met with several laser cooling techniques, as for ex- 
ample Sisyphus cooling HO] , velocity selective coher- 
ent population trapping 111, Raman cooling [12], Raman 
sideband cooling |13| , or Doppler cooling on narrow lines 
|14j . The second challenge is the implementation of an 
efficient trapping scheme that allows for accumulation 
of atoms at high density in a particular region |15| 116] . 
The third, and most severe challenge is to avoid the detri- 
mental effects of the laser cooling photons, which impede 
the required density increase. One such effect are light- 



assisted inelastic collisions, which lead to loss [T71 118). 
Another is the reabsorption of photons scattered dur- 
ing laser cooling [TS], which leads to an effective repul- 
sion between the atoms and to heating of atoms in the 
lowest energy states. Both effects increase with density 
and make it impossible to reach quantum degeneracy. 
For low phase-space density samples, this challenge has 
been overcome by rendering the atoms transparent to 
laser cooling photons 20- 22J or by decreasing the pho- 
ton scattering rate below the frequency of a confining 
trap [3 [33] . It has also been proposed to reduce re- 
absorption by dimensional reduction of the sample [7]- 
The solutions to the three challenges implemented so far 
are insufficient to reach quantum degeneracy. The high- 
est phase-space densities ever attained are one order of 
magnitude too low [HI [24] . Surprisingly, this last order 
of magnitude has been an insurmountable obstacle for a 
decade. 

In this Letter, we present an experiment hat overcomes 
all three challenges and creates a BEC of strontium by 
laser cooling. Our scheme is based on the combina- 
tion of three techniques, favored by the properties of 
this element, and does not rely on evaporation. Stron- 
tium possesses a transition with such a narrow linewidth 
(r/2-7r = 7.4 kHz) that simple Doppler cooling can reach 
temperatures down to 350 nK [TU [2SJ [53] ■ Using this 
transition, we prepare a laser cooled sample of 10 7 84 Sr 
atoms in a large "reservoir" dipole trap. To avoid the 
detrimental effects of laser cooling photons, we render 
atoms transparent for these photons in a small spatial 
region within the laser cooled cloud. Transparency is in- 
duced by a light shift on the optically excited state of the 
laser cooling transition. In the region of transparency, 
we are able to increase the density of the gas, by accu- 
mulating atoms in an additional, small "dimple" dipole 
trap [HI [16] . Atoms in the dimple thermalize with the 
reservoir of laser-cooled atoms by elastic collisions and 
form a BEC. 

The details of our scheme are shown in Fig. [T] Based 
on our previous work [2"TH2"9"] . we use several stages of 
laser cooling to prepare a sample of 84 Sr atoms in the 



2 



(a) 



I 



transparency" 
beam 



(b) 




688 nm 
4.4 MHz 



689 nm 
7.4 kHz 



cooling light 



'transparency" 
transition 



cooling 
transition 



y <8> > 



light shift by 
"transparency" beam 





Figure 1: Scheme to reach quantum degeneracy by laser cool- 
ing, (a) A cloud of atoms is confined in a deep reservoir dipole 
trap and exposed to a single laser cooling beam (red arrow). 
Atoms are rendered transparent by a "transparency" laser 
beam (green arrow) and accumulate in a dimple dipole trap 
by elastic collisions, (b) Level scheme showing the laser cool- 
ing transition and the transparency transition, (c) Potential 
experienced by 1 So ground-state atoms and atoms excited to 
the 3 Pi state. The transparency laser induces a light shift 
on the 3 Pi state, which tunes the atoms out of resonance 
with laser cooling photons, (d) to (f) Absorption images of 
the atomic cloud recorded using the laser cooling transition. 
The images show the cloud from above and demonstrate the 
effect of the transparency laser (e) and the dimple (f). (d) is 
a reference image without these two laser beams. 



reservoir trap [30] . The trap consists of an infrared 
laser beam (wavelength 1065 nm) propagating horizon- 
tally (^-direction). The beam profile is strongly elliptic, 
with a beam waist of 300 /im in the horizontal direc- 
tion (y-direction) and 17 ^m along the field of gravity 
(z-direction). The depth of the reservoir trap is kept 
constant at ks x 9/iK. After preparation of the sam- 
ple, another laser cooling stage is performed on the nar- 
row Sq — 3 Pi intercombination line, using a single laser 
beam propagating vertically upwards. The detuning of 
the laser cooling beam from resonance is ~ —2 T and the 
peak intensity is 0.15 /j,W/cm 2 , which is 0.05 of the tran- 
sition's saturation intensity. These parameters result in 
a photon scattering rate of ~ 70s -1 . At this point, the 
ultracold gas contains 9 x 10 6 atoms at a temperature of 
900 nK. 

To render the atoms transparent to cooling light in a 
central region of the laser cooled cloud, we induce a light 



shift on the 3 Pi state, using a "transparency" laser beam 
15 GHz blue detuned to the 3 P X - 3 Si transition [30] . 
This beam propagates downwards under a small angle 
of 15° to vertical and has a beam waist of 26 /.tm in the 
plane of the reservoir trap (xy-plane) . The beam has a 
peak intensity of 0.7kW/cm 2 . It upshifts the 3 Pi state 
by more than 10 MHz and also influences the nearest 
molecular level tied to the 3 Pi state significantly |5Ull31j . 
Related schemes of light-shift engineering were used to 
image the density distribution of atoms |32| I33j , to im- 
prove spectroscopy [31] , or to enhance loading of dipole 
traps [2D 122] ■ To demonstrate the effect of the trans- 
parency laser beam, we take absorption images of the 
cloud on the laser cooling transition. Figure [ljd) shows 
a reference image without the transparency beam. In 
presence of this laser beam, atoms in the central part of 
the cloud are transparent for the probe beam, as can be 
seen in Fig. [TJe) . 

To increase the density of the cloud, a dimple trap 
is added to the system. It consists of an infrared laser 
beam (wavelength 1065 nm) propagating upwards under 
a small angle of 22° to vertical and crossing the laser 
cooled cloud in the region of transparency. In the plane 
of the reservoir trap, the dimple beam has a waist of 
22 /im. The dimple is ramped to a depth of ks x 2.6 //K, 
where it has trap oscillation frequencies of 250 Hz in the 
horizontal plane. Confinement in the vertical direction 
is only provided by the reservoir trap and results in a 
vertical trap oscillation frequency of 600 Hz. Figure [ljf) 
shows a demonstration of the dimple trap in absence of 
the transparency beam. The density in the region of the 
dimple increases substantially. However, with the dimple 
alone no BEC is formed because of photon reabsorption. 

The combination of the transparency laser beam and 
the dimple trap leads to BEC. Starting from the laser 
cooled cloud held in the reservoir trap, we switch on the 
transparency laser beam and ramp the dimple trap to a 
depth of ks x 2.6 [iK. The potentials of the 1 Sq and 3 Pi 
states in this situation are shown in Fig.[T|c). Atoms ac- 
cumulate in the dimple without being disturbed by pho- 
ton scattering. Elastic collisions thermalize atoms in the 
dimple with the laser cooled reservoir. The phase-space 
density in the dimple increases and a BEC emerges. 

We detect the BEC by taking absorption images 24 ms 
after switching off all laser beams. Figure [2]ja) shows 
the momentum distribution 20 ms after switching on the 
transparency beam, which is well described by a ther- 
mal distribution. By contrast, we observe that 140 ms 
later, an additional, central elliptical feature has devel- 
oped; see Fig. [2^b) . This is the hallmark of the BEC. 
Although clearly present, the BEC is not very well visi- 
ble in Fig.^b), because it is shrouded by 8 x 10 6 thermal 
atoms originating from the reservoir. To show the BEC 
with higher contrast, we have developed a background 
reduction technique. We remove the reservoir atoms by 
an intense flash of light on the 1 Sq — 3 Pi transition ap- 
plied for 10 ms. Atoms in the region of transparency 
remain unaffected by this flash. Only 5 x 10 5 thermal 
atoms in the dimple remain and the BEC stands out 
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Figure 2: Creation of a BEC by laser cooling. Shown are time-of-fiight absorption images and integrated density profiles of 
the atomic cloud for different times t after the transparency laser has been switched on, recorded after 24 ms of free expansion, 
(a) and (b) The appearance of an elliptic core at t — 160ms indicates the creation of a BEC. (c) Same as in (b), but to 
increase the visibility of the BEC, atoms in the reservoir trap were removed before the image was taken. The fits (blue lines) 
consist of Gaussian distributions to describe the thermal background and an integrated Thomas-Fermi distribution describing 
the BEC. The red lines show the component of the fit corresponding to the thermal background. The x'y'-plane is rotated by 
45° around the z-axis with respect to the xy-plane and the field of view of the absorption images is 2 mm x 1.4 mm. 



clearly; see Fig. |2jc) . We use this background reduction 
technique only for demonstration purposes, but not for 
measuring atom numbers or temperatures. 

Quantitative data on our experiment are obtained by 
two-dimensional fits to time-of-flight absorption images 
[3D]. The atom number of the thermal cloud and of the 
BEC are extracted from fits to 24- ms expansion images, 
consisting of Gaussian distributions describing the ther- 
mal background and an integrated Thomas-Fermi distri- 
bution describing the BEC. Further absorption images 
taken after 4 ms expansion time are used to determine 
atom number and temperature of the gas in the reservoir 
and the dimple, respectively. 

We now analyze the dynamics of the system after the 
transparency laser beam has been switched on. As we 
increase the dimple strength to its final depth in 10 ms, 
10 6 atoms accumulate in it and the temperature of the 
dimple gas increases; see Figs. [3|a) and (b). During the 
next ~ 100 ms the dimple gas thermalizes with the reser- 
voir gas by elastic collisions [351 ES] ■ The temperature 
of the reservoir gas is hereby not increased, since the 
energy transferred to it is dissipated by laser cooling. 
We carefully check that evaporation is negligible even 
for the highest temperatures of the gas [30] . Already 
after 60 ms a BEC is detected. Its atom number satu- 
rates at 1.1 x 10 5 after 150ms, as shown in Fig. |3|c). 
The atom number in the reservoir decreases slightly, ini- 
tially because of migration into the dimple and on longer 
timescales because of light assisted loss processes in the 
laser cooled cloud. 

To demonstrate the power of our laser cooling scheme, 
we repeatedly destroy the BEC and let it reform (Fig. |4| . 
To destroy the BEC, we pulse the dimple trap depth to 



kg x 15 /iK for 2 ms, which increases the temperature of 
the dimple gas by a factor two. We follow the evolu- 
tion of the BEC atom number while the heating pulse is 
applied every 200 ms. A new BEC starts forming a few 
10 ms after each heating pulse for more than 30 pulses. 
We find that the observed decrease in the BEC size from 
pulse to pulse stems from the reduction of the total atom 
number in the system. 

To clarify the role laser cooling plays in our scheme, we 
perform a variation of the experiment. Here, we switch 
off the laser cooling beam before ramping up the dim- 
ple and we do not use the transparency beam. Heat 
released while ramping up the dimple or after a heating 
pulse is again distributed from the dimple to the whole 
system by elastic collisions, but this time not dissipated 
by laser cooling. Since the reservoir gas has a ten times 
higher atom number than the dimple gas, the temper- 
ature after thermalization is only increased by a small 
amount. If the final temperature in the dimple is below 
the critical temperature, a BEC is formed. This scheme 
resembles the formation of a BEC by trap deformation, 
as demonstrated in [15] using a sample of atoms cooled 
by evaporation. We test the performance of this BEC 
creation scheme again by repeated heating pulses. We 
can detect a BEC after at most five pulses. For more 
pulses, the temperature of the gas in the dimple remains 
too high to allow the formation of a BEC. This poor be- 
havior stands in stark contrast to the resilience of BEC 
formation to heating, if the system is continuously laser 
cooled. 

The ability to reach the quantum degenerate regime 
by laser cooling has many exciting prospects. This 
method can be applied to any element possessing a laser 
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Figure 3: Characterization of the BEC formation process 
after the transparency laser is switched on. Shown are the 
evolution of the atom number in the dimple and the reservoir 
(a), the evolution of temperature in these regions (b) and 
the BEC atom number (c). During the first 10 ms of this 
evolution, the dimple trap is ramped on. After 60 ms a BEC 
is detected. 
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Figure 4: Repeated destruction and reformation of the BEC. 
Shown is the evolution of BEC atom number while the BEC 
is destroyed every 200 ms (arrows) by suddenly increasing the 
depth of the dimple trap. If the system is laser cooled, the 
BEC atom number quickly increases again, which is shown 
here for up to 18 cycles (red circles). Without laser cooling, 
a BEC is detectable for at most five heating cycles, of which 
the first three are shown here (blue squares). 
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cooling transition with a linewidth in the kHz range 
and suitable collision properties. Besides strontium this 
encompasses several lanthanides (37J ISH] • The technique 
can also cool fermions to quantum degeneracy and it 
can be extended to sympathetic cooling in mixtures 
of isotopes or elements. Another tantalizing prospect 
enabled by variations of our techniques, is the realization 
of a continuous atom laser, which converts a thermal 
beam into a laser-like beam of atoms. To realize such a 
device, our scheme needs to be extended in two ways. 
First, the reservoir of laser cooled atoms needs to be 
replenished, for example by sending a thermal beam 
of atoms onto a part of the reservoir with sufficiently 
high cooling laser intensity to allow capture of these 
atoms. Second, a continuous beam of condensed atoms 
needs to be outcoupled. Using magnetic species such 
as dysprosium or erbium, outcoupling from the BEC 
is possible by changing the internal state and thereby 
the magnetic force on the atoms [2H1 SO] • Alternatively, 
the reservoir can be connected to an outcoupling dipolc 
trap, creating a narrow channel where atoms can escape 
without further interaction with the cooling light [41]. 
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Supplemental Material 

This Supplemental Material contains in-depth infor- 
mation about our approach to reach quantum degener- 



acy by laser cooling. In Sec. 
perimental sequence. In Sec. 



I] we give details on the ex- 
Tllwe discuss the parameter 



dependence of BEC creation and properties of the BEC. 
In Sec. |III| we show that evaporative cooling does not 
play a role in reaching quantum degeneracy. In Sec. |IV| 
we analyze the shifting of atomic and molecular transi- 
tions by the transparency beam. In Sec. [V] we model the 
density distributions of the BEC and the thermal cloud. 



I. EXPERIMENTAL SEQUENCE 

In the following, we discuss in detail the experimental 
sequence with which we obtain BEC by laser cooling. We 
first describe the preparation of a sample of ultracold 

Sr atoms in the reservoir dipole trap, which follows 
closely the procedure used in our previous work |27H29] . 
Then we give details on the additional steps we take to 
produce a BEC. These details include information on 
the transparency beam, the dimple dipole trap and the 
conditions of the final laser cooling stage. We end with a 
discussion of our data acquisition and analysis method. 



A. Sample preparation 

Isotope choice — Of the three bosonic strontium 
isotopes, Sr is best suited for our experiment, since its 
scattering length of as4 = 124 ao (with ao the Bohr ra- 
dius) allows efficient thermalization by elastic collisions. 
The other bosonic isotopes are unsuitable for our task. 
88 Sr has a negligible scattering length of ass = — 1.4 do 
and therefore does not thermalize. 86 Sr suffers from 
three-body inelastic loss because of a very large scat- 
tering length ase = 800 ao- Unfortunately 84 Sr has a low 
natural abundance of only 0.56%. To obtain a large sam- 
ple, we accumulate metastable state atoms in a magnetic 
trap, as described below. 

Atomic beam and blue MOT — A sample of 
strontium metal with natural isotopic composition is 
heated in an oven under vacuum to about 600 °C in 
order to sublimate strontium atoms. The atoms form 
an atomic beam after escaping the oven through a bun- 
dle of microtubes. 84 Sr atoms in the atomic beam are 
transversally cooled, Zeeman-slowed, and captured in 
a "blue" magneto-optical trap (MOT) with laser light 
red-detuned to the 1 S'o — 1 P\ transition at 461 nm; see 

2 /2tt = 



Fig. SI This transition has a linewidth of Tbi u 
30.5 MHz, corresponding to a Doppler temperature of 
T DM uc = hr Uuc /(2k B ) = 720 fiK. The MOT uses a 
quadrupole magnetic field with vertically oriented axis 
and a vertical field gradient of 55G/cm. The MOT 
beams have waists of 5 mm and peak intensities of 
10mW/cm 2 , which is a quarter of the saturation inten- 



sity I S at,biuc = 7r/icr b iuc/3Ag luc = 40.7mW/cm 2 . The 
detuning of the MOT light to resonance is —30 MHz. 

Metastable state reservoir — The blue MOT cy- 
cle is not completely closed, as atoms in the excited 1 P\ 
state can decay into the metastable and magnetic 3 P 2 
state with a branching ratio of 1:150 000. Atoms in this 
state can be trapped in the quadrupole magnetic field 
of the MOT. Since the lifetime of magnetically trapped 
3 P2-state atoms of ~ 30 s is about three orders of mag- 
nitude longer than the leak time scale of the blue MOT, 
metastable state atoms accumulate in the magnetic trap. 
We operate the MOT until about 10 s atoms are accu- 
mulated in this "metastable state reservoir" , which takes 
~ 10 s. This accumulation of atoms allows us to use the 

Sr isotope, despite its low natural abundance. 

Red MOT — To increase the phase-space density 
of the sample, we use a narrow-line "red" MOT op- 
erated on the 1 iSo — 3 Pi intercombination line, which 
has a wavelength of A rc d = 689 nm and a linewidth of 
r rod /27r = 7.4 kHz. Our red MOT laser has a linewidth 
of about 2 kHz and an absolute stability better than 
500 Hz. The horizontal (upward, downward propagat- 
ing) red MOT beams have a waist of 3.1mm (2.9 mm, 
3.1mm). The Doppler temperature of the red MOT 
Po.red = 180 nK is comparable to the recoil tempera- 
ture T r = h 2 k 2 /(k B m) = 460 nK, where k = 2ir/\ rcd is 
the wave vector of the cooling light, and m is the atom's 
mass. The temperature can approach T r /2 when reduc- 
ing the cooling light intensity to the saturation inten- 
sity J S at, red = 3 /jW/cm 2 [35]. Not only the temperature 
but also the atom number decreases, when decreasing 
the cooling light intensity, and a compromise between 
atom number and temperature has to be chosen. We 
typically reach temperatures of 800 nK with samples of 
about 10 7 atoms. 
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Figure SI: Level scheme of strontium with all relevant states 
and transitions. The dotted arrows show the decay path of 
atoms from the 1 Pi state into the 3 P2 state. 
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To load the red MOT, metastable state atoms are 
transferred into the 1 S'o ground state by optical pump- 
ing via the 5s5d 3 D 2 state. The temperature of the sam- 
ple recovered from the metastable state reservoir is on 
the order of the Doppler temperature of the blue MOT, 
Tb,biue = 720 /iK, which is three orders of magnitude 
higher than the final temperature reached with the red 
MOT. We initially use red MOT parameters that allow 
to capture such a high temperature sample and then 
ramp these parameters to conditions in which a high 
phase space density is reached. The intensity and fre- 
quency of the MOT light during the capture phase and 
the subsequent ramps are given in Figs. |S2|^a) and (b). 
During the capture phase, the red MOT laser beams are 
frequency-modulated to increase the capture velocity of 
the MOT. In addition, the quadrupole magnetic field 
gradient is suddenly lowered to 1.15(5) G/cm in the ver- 
tical direction, which increases the capture volume, and 
then maintained at this value for the remainder of the 
experimental sequence. The detuning of the MOT light 
is reduced to compress the atomic cloud, and the inten- 
sity is reduced to lower the temperature. This provides 
optimal conditions for loading the reservoir dipole trap. 

Reservoir dipole trap — The reservoir dipole trap 
is a large-volume dipole trap consisting of a horizon- 
tally propagating beam. This beam is derived from a 
5-W fiber laser operating at a wavelength of 1065 nm 
with a linewidth of 0.5 nm (IPG ytterbium fiber laser, 
model YLD-5-1064-LP). The beam contains up to 2 W 
of power and is linearly polarized in the vertical direc- 
tion. Cylindrical lenses are used to create a horizon- 
tally elongated beam profile. The beam has a vertical 
(horizontal) waist of 17.3(2) /im (298(3) /mi), yielding 
an aspect ratio of about 1:17. The trap frequencies are 
600(6) Hz (35.3(4) Hz, 5.78(6) Hz) in the vertical (trans- 
verse horizontal, axial horizontal) direction at a power 
of 1760(35) mW as used in the experiment. Including 
gravitational sagging, the potential depth in the vertical 
direction is ~ kg x 9 /iK. The errors given here are domi- 
nated by the uncertainty in the measurement of the laser 
beam powers, which we assume to be 2%. The uncer- 
tainty in the trap oscillation frequencies is much smaller, 
and other effects, such as transmission losses of the glass 
cell and gravitational sagging in the dipole trap, are ac- 
counted for. 

We increase the axial frequency slightly by a sec- 
ond, near-vertical beam with 297(3) /im waist and cir- 
cular polarization, using a second fiber laser of the same 
model as mentioned above. At the power of 740(15) mW 
used in the experiment, it has a trap depth of only 
k B x 0.29(1) fiK and increases the trap frequency in the 
x-direction to 8.14(8) Hz. This increased axial trap fre- 
quency slightly increases the density of the sample and 
reduces the timescale of BEC formation. It has no fur- 
ther effects, and otherwise identical results are obtained 
without this additional beam. 

Loading of the reservoir dipole trap — When 
the intensity of the MOT beams is low, gravity plays 
an important role in the MOT dynamics. The atoms 



are pulled towards the lower part of an ellipsoid of equal 
magnetic field magnitude located around the quadrupole 
magnetic field center. On this ellipsoid, the atoms are 
in resonance with the cooling light and levitated by it, 
giving the MOT a pancake shape, roughly matched by 
the shape of the reservoir dipole trap. The dipole trap 
is located 600(100) /jm below the center of quadrupole 
field, where the B-field has a magnitude of 75(10) mG 
and is oriented nearly vertically. To overlap the MOT 
with the dipole trap, we carefully adjust the MOT laser 
detuning. Figures S2 'c) and (d) show the temperature 
and atom number evolution of the atomic cloud, while 
it is loaded into the reservoir dipole trap. At the end- 
point of the ramp, the intensities of the horizontal (up- 
ward, downward propagating) MOT beams are 0.6 7 sat 
(0.5 / sa t, / S at)- The detuning from the cr + -cooling tran- 
sition is about —20 kHz. With the given non-zero B-field 
at the position of the atoms, this corresponds to a detun- 
ing of about 150 kHz from the unperturbed 7r-transition. 

Laser cooling of the atoms in the dipole trap is influ- 
enced by Zeeman and light shifts. The vertical gradient 
of 1.15 G/cm translates to a negligible shift of 600 /iG 
across the sample, assuming an estimated vertical size of 
5 /im. The horizontal gradient is 0.63 G/cm, and the hor- 
izontal radius of the cloud is rather large: about 200 /im, 
giving a shift of 12 mG between the center and the out- 
sides, which corresponds to a significant frequency shift 
of 25 kHz. The dipole trap induces a non- negligible dif- 
ferential AC Stark shift on the 1 So and 3 Pi states. We 
can reduce this shift by choosing the optimized polariza- 
tions for the horizontal and vertical dipole trap beams 
mentioned above, such that the magnitude of the shift 
is about 10 kHz [26]. Assuming that the atoms explore 
one tenth of the trap depth, we obtain a shift across the 
sample of about 1kHz, which is negligible. 

About 9 x 10 6 atoms are captured in the reservoir 
dipole trap at a temperature of 830 nK. The achievable 
temperature is density-dependent. A reduction of the 
atom number by a factor two leads to a temperature re- 
duction of ~ 100 nK. Over the course of one hour, the 
variation of atom number between experimental cycles is 
about 1%, and the temperature variation is below 10 nK. 
After loading of the dipole trap, the intensity of the up- 
ward cooling beam is reduced by a factor of 10, and 
the horizontal beams are turned off. Such a sample is 
the starting point for the subsequent laser cooling into 
quantum degeneracy. 



B. Transparency beam 

Immediately after completion of the dipole trap load- 
ing, the transparency beam is turned on. This beam is 
derived from a free-running master diode laser with a fre- 
quency stability of order 100 MHz/day. The frequency is 
blue detuned by 15 GHz from the 3 Pi — 3 Si transition, 
and constantly monitored by a wavemeter. The beam is 
circularly polarized and focused onto the center of the 
atomic cloud with a waist of 26.2(3) /im, again calcu- 
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Figure S2: Narrow-line MOT and dipole trap loading. The narrow-line cooling consists of five consecutive stages, labeled 
MOT 1 through V: a 50-ms capture MOT, during which the repumping of atoms from the Z P2 to the 1 So state takes place, 
a compression MOT of 100 ms, two cooling MOT stages of 100 ms and 250 ms, and a 100-ms hold time at final parameters. 
The end of the MOT phase is chosen as the origin of the time axis, (a) Intensity of the MOT beams in units of the saturation 
intensity J sa t = 3/iW/cm 2 . The intensity of the horizontal, upward, and downward MOT beams can be set independently. 
The downward MOT beam has little effect and is turned off before the other ones (see inset), (b) Frequency of the MOT 
beams, given as detuning from the Zeeman-shifted a + transition. The light is frequency-modulated in the first two MOT 
stages. For clarity, only every tenth line of the resulting frequency comb is shown here. Atom number (c) and temperature 
(d) drop significantly during the end of MOT stage IV. The drop in atom number happens as the restoring force of the MOT 
becomes too weak to support atoms against gravity, and all atoms outside the dipole trap are lost. The last 100 ms of hold 
time do not reduce the temperature further, but this stage is crucial to increase the density of atoms in the dipole trap. 



lated from trap frequency measurements. The power 
used is 7.5(2) mW, translating to a peak intensity of 
7 x 10 5 mW/cm 2 = 7MW/m 2 . The beam propagates 
downwards, at an angle of 15° to vertical due to geomet- 
rical restrictions. 

The transparency beam has a strong influence on the 
3 Pi state. The differential AC Stark shift of the 1 S - 3 Pi 
cooling transition is on the order of +10 MHz, such that 
the cooling light is red-detuned by > 1000 r ro d for atoms 



at the center of the transparency beam. In this way, 
atoms illuminated by the transparency beam are trans- 
parent to laser cooling photons. Note that this scheme 
is drastically different to a scenario in which the cool- 
ing beam would contain a small dark spot imaged onto 
the dimple region: in this case, atoms in the dimple re- 
gion could still absorb cooling light scattered by atoms 
in the reservoir. Our method differs also from the dark 
spot MOT technique, since we do not change the inter- 
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Table SI: Parameters of the dipole trap beams used for the experiments described in the Letter. Gravitational sagging is taken 
into account in the calculation of the horizontal dipole trap depth. 



beam 


waist x 


waist y 


waist z 
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U/k B 


u 


fv 


/. 




W 


(lim) 


W 


(mW) 


(/iK) 


(Hz) 


(Hz) 


(Hz) 


horizontal 




298(3) 


17.3(2) 


1760(35) 


9.2(2) 


5.78(6) 


35.3(4) 


600(6) 


vertical 


297(3) 


297(3) 




740(15) 


0.29(1) 


5.73(6) 


5.73(6) 


~ 


dimple 


22.4(2) 


22.4(2) 




38.3(8) 


2.6(1) 


228(2) 


228(2) 


~ 



nal state of the atom in the region of transparency. A 
careful analysis of the magnitude of the AC Stark shift 



of the 3 Pi state is given in Sec. IV 



The transparency beam also creates an attractive 
trapping potential for the x 5o state with a depth of 
ks x 0.5 /iK. Because of similar beam orientation and 
waists, the potential created by the transparency beam 
resembles the potential created by the dimple beam, but 
has only about 20% of its depth. 



C. Cooling light 

The cooling light consists of an upward propagating, 
circularly polarized beam, red detuned by about 15 kHz 
from the a + 1 So— 3 Pi transition. It has a power of 20 nW 
and a peak intensity of 0.15/iW/cm 2 , corresponding to 
0.05 I S at- The light has the same source as the light used 
for the narrow-line MOT. Any desired cooling rate of up 
to many 100 nK/ms can be achieved by a suitable combi- 
nation of detuning and intensity, where a larger detuning 
can be compensated by an increased intensity. We find 
a fixed relation between cooling rate and induced decay 
rate, which is independent on the combination of detun- 
ing and intensity. There is a lower temperature limit 
of about 450 nK, below which cooling is accompanied 
by rapidly increasing loss rates. For the experiments 
described in the Letter, we use only one upward prop- 
agating cooling beam, and we find that the addition of 
further cooling beams from other directions does not im- 
prove the performance. 



D. Dimple 

The local increase in density is facilitated by the dim- 
ple beam. This beam is aligned almost vertically, with 
an angle of 22° (37°) towards vertical (the transparency 
beam) due to geometrical restrictions. It is derived from 
the same laser source as the vertical dipole trap, has cir- 
cular polarization and a waist of 22.4(2) /ini. The dim- 
ple beam is centered in the plane of the reservoir trap 
with the transparency beam to within 5 /im. At a power 
of 38.3(8) mW, the dimple provides horizontal trap fre- 
quencies of 228(2) Hz and has negligible trap frequencies 
in the vertical direction. In presence of the dimple, we 
refer only to the region outside of the dimple as the reser- 
voir. 

The dimple trap is set to a small depth of fcg x 



0.15 /iK at the beginning of the experimental sequence 
and ramped in 10 ms to a depth of ks x 2.6 /iK after 
the transparency beam is switched on. The ramp speed 
is not adiabatic with respect to the horizontal trap fre- 
quencies of the reservoir trap, and atoms from the reser- 
voir continue to accumulate in the dimple for 10 ms after 
the ramp. At this point, the temperature of the gas in 
the dimple is twice the temperature of the reservoir gas. 
Thermal contact with the laser-cooled reservoir lowers 
the temperature on a timescale of ~ 100 ms. In thermal 
equilibrium, the dimple leads to a peak density increase 
by a factor of 30 compared to the reservoir. This density 
increase is the origin of the gain in phase-space density. 

The BEC phase transition is observed 50 ms after 
ramping the dimple to high power. The BEC grows 
to slightly more than 10 5 atoms after another 100 ms. 
The local harmonic potential confining the BEC is dom- 
inated horizontally by the dimple trap and vertically 
by the horizontal dipole trap, leading to trap frequen- 
cies of 228(2) Hz in the horizontal plane and 600(6) Hz 
in the vertical direction. The BEC is pancake-shaped 
with a horizontal (vertical) Thomas- Fermi radius of 5 /im 
(2/im). 



E. Absorption imaging and data analysis 

We use time-of-flight absorption images to deduce all 
relevant information from our atomic samples. We are 
interested in the following quantities: atom number in 
the dimple and in the reservoir, temperature of the dim- 
ple and of the reservoir, and number of atoms in the 
condensate. 

The temperature and atom number of the dimple 
and reservoir regions are measured in absorption images 
taken after 4 ms of free expansion; see Fig. |S3[a). After 
this short time, the clouds from the two regions can still 
be clearly distinguished in the horizontal direction (in- 
direction), but have expanded well beyond their initial 
size in the vertical direction (^-direction) to allow ther- 
mometry. Since the atomic clouds are very dense and 
nearly opaque to resonant imaging light, we image at a 
detuning of 48 MHz 1.5 T, which reduces the absorp- 
tion cross section by a factor of about 12. We employ a 
2D double-Gaussian fit to the data. The temperature is 
derived from the vertical expansion only. 

The overall atom number and the BEC atom number 
are determined from absorption images taken after a free 
expansion time of 24ms; see Fig. S3 b). We fit the data 
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x' (mm) x' (mm) 

Figure S3: Typical absorption images used to determine tem- 
peratures and atom numbers, (a) Absorption image after a 
short expansion time t cx = 4 ms. A two-dimensional double- 
Gaussian fit is used to extract atom number and tempera- 
ture individually for dimple and reservoir, (b) Absorption 
image after a long expansion time i C x = 24 ms. A tri-modal 
two-dimensional fit, consisting of two Gaussians and an in- 
tegrated Thomas-Fermi distribution, are used to extract the 
atom number of the thermal component and the BEC. The 
panels below the absorption images show the density distri- 
bution integrated along the ^-direction (black) and the fits 
(blue) . The red curves show the contributions of the reservoir 
(a) and of the thermal component (b). 



dependence of the total atom number N, at these low- 
est achievable temperatures. We find that iVeEC oc N. 
This behavior can be explained by the model described 
in Sec. El 

Trap depth — We independently vary the trap 
depth of horizontal beam and dimple around the val- 
ues presented in the Letter. A variation in trap depth 
of the horizontal beam by a factor of 1.7 changes the 
BEC atom number by at most 15%, and a variation of 
dimple power by a factor of 2.5 changes the atom num- 
ber by at most 30%. A broad, global maximum of the 
BEC atom number exists, which is where we perform 
our experiments. 

An important constraint on the trap depth is the de- 
sire to avoid evaporation. This constraint is fulfilled by 
the conditions chosen; see Sec. |III| for details. 

Dimple — We find that the BEC atom number de- 
pends only weakly on the dimple size. We vary the waist 
of the dimple beam between 20 and 50 /im, but find no 
appreciable change in BEC atom number. The BEC for- 
mation time, however, is increased for larger waists. We 
also vary the ramp-up time of the dimple between and 
100 ms, but find no influence on temperature or BEC 
atom number after a 250-ms equilibration time. 



B. Transparency beam 



with a 2D distribution consisting of two Gaussians for 
the thermal atoms and an integrated Thomas-Fermi dis- 
tributions for the BEC. Two Gaussians are used to take 
into account the two different sources of thermal atoms, 
the reservoir and the dimple. If only one Gaussian is 
employed to describe all thermal atoms, the BEC atom 
number deduced from the Thomas-Fermi part of the fit 
increases by about 50%. 

All measurements are performed three times, and the 
average values and statistical error bars of atom number 
and temperature are given in Figs. 3 and 4 of the Letter. 



II. CHARACTERIZATION OF LASER-COOLED 
BECS 

The BEC formation depends on various parameters. 
In this section, we will explore the influence of these 
parameters on the system and show that our method 
of laser cooling to quantum degeneracy works in a very 
broad range of parameters. 



A. BEC atom number 

Total atom number and temperature — The 

lowest achievable temperature depends on the total atom 
number in the system as a consequence of red MOT dy- 
namics. We measure the BEC atom number ATbec m 



The transparency beam is the key novelty of our work. 
A detailed study of the non-trivial shifting of the laser 



cooling transition can be found in Sec. IV 

Frequency — The transparency beam is employed 
to locally shift the energy of the 3 Pi state by about 
+10 MHz. This would in principle be possible with light 
blue-detuned to any transition originating from the 3 Pi 
state. We chose to work in the vicinity of the 3 Pi — 3 Si 
transition at 688 nm for two reasons: The availability 
of diode lasers, and the negligible influence on the 1 Sq 
state. 

In our experiment, the transparency beam is blue- 
detuned by about 15 GHz from the 3 Pi — 3 Si transition. 
In a series of measurements, we set the detuning to dif- 
ferent values and vary the intensity of the transparency 
beam while searching for the appearance of a BEC. We 
never observe BEC formation for frequencies around or 
red-detuned to the 3 Pi — 3 Si transition. We observe 
BEC formation only in a frequency range between 6 
and 30 GHz blue detuning, where the upper bound of 
30 GHz is probably limited by the available laser power 
of 10 mW. 

Intensity — We observe the formation of a BEC 
only above a certain critical intensity I c of the trans- 
parency beam, which depends on its frequency. The 
BEC atom number quickly grows for larger intensities, 
and saturates at about 2I C . The experiment is per- 
formed at around 10 I c , corresponding to a peak intensity 
of 0.7kW/cm 2 . 

Spectral filtering — The 3 Pi — 3 Si transition is 
1.428 nm (or 902 GHz or 30 cm" 1 ) away from the 1 S ~ 
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3 Pi intercombination transition, and off-resonant excita- 
tion of the intercombination transition is negligible. The 
light originating from the laser diode, however, contains 
incoherent fluorescence light, also known as residual am- 
plified spontaneous emission (ASE). This light covers a 
broad spectrum with a width of about 20 THz, which in- 
cludes the 1 S'o — 3 Pi transition. A spectrometer is used 
to estimate the spectral power of the incoherent back- 
ground Pase to be on the order of 10~ 12 Po in a 1kHz 
wide band 1 THz away from the laser line, where Pq is the 
power in the laser line. The scattering of these incoher- 
ent photons limits the lifetime of a BEC to 800 ms when 
illuminated by unfiltered transparency light. Spectral 
filtering of the light allows us to reduce the amount of 
resonant light by a factor of 500, leading to an increased 
lifetime of 10 s. The lifetime measurement is performed 
with a pure BEC obtained by standard evaporation, and 
its lifetime is limited to 10 s by inelastic collisions. 

Waist — The transparency beam needs to be well- 
aligned with the dimple beam, and it needs to cover the 
BEC entirely. We vary the waist of the transparency 
beam between 10 and 55 /xm. For comparison, the dimple 
waist is 22 /jm and the Thomas-Fermi radius of the BEC 
in the xy-plane is 5 /xm. Given an appropriate adjust- 
ment of the intensity, we can create laser-cooled BECs 
within the entire range of transparency beam waists ex- 
amined. The purification method described in Fig. 2(c), 
during which reservoir atoms are selectively removed, be- 
comes increasingly inefficient for larger beam sizes, as 
also parts of the reservoir are shielded. A larger beam 
size makes the system insensitive to misalignments of the 
beams, which we verify experimentally. 

Local density increase — One might speculate that 
the transparency beam leads to an increase in atom den- 
sity purely by elimination of photon emission and re- 
absorption cycles, which act as an effective repulsion. 
This is not the case. We do observe a small density in- 
crease induced by the transparency beam, but it can be 
explained entirely by its trapping potential for atoms in 
the 1 Sq state. 



C. Timescales 

Collision rate and thermalization — The peak 
scattering rate of thermal atoms in the dimple region is 
3100 s" 1 , found at the edge of the BEC (see Sec|vj). The 
scattering rate decreases with distance from the center of 
the dimple, down to 50s _1 outside of the dimple region. 
We observe that thermalization happens on a timescale 
of about 100 ms, which corresponds to about 5 collisions; 
compare Fig. 3(b). This number matches well with the 
expectation of about 3 collisions required for thermaliza- 
tion. 

BEC formation time — We study the time required 
for the BEC to form after the dimple has been ramped 
up. For a BEC to be created, the atom number in the 
dimple has to be high enough and the temperature below 
a critical temperature for the given atom number. We 



find that the timescale required to accumulate atoms in 
the dimple is much shorter than the timescale required to 
thermalizc the dimple gas with the reservoir; see Fig. 3. 
Thus, the BEC formation time is limited by the thermal- 
ization time scale. This time scale is linked to the elastic 
collision rate, which depends on the density of the sam- 
ple. To test this relation, we vary the density by two 
means, changing the reservoir atom number and chang- 
ing the reservoir trap frequency. We observe an increase 
in formation time as we decrease the atom number, down 
to a minimum atom number of about 1 x 10 6 atoms, be- 
low which no BEC forms in our trap geometry. Using 
the trap parameters of the experiment described in the 
Letter, the fastest formation time is about 50 ms for the 
maximum achievable atom number of 10 x 10 6 . An in- 
crease of the axial frequency from 8 Hz to 40 Hz reduces 
the formation time to about 20 ms. The timescale of 
BEC formation (following abrupt changes of the ther- 
mal distribution) was studied in previous experiments 
[351 136] and found to depend on the scattering rate, in 
agreement with our measurements. 

Reservoir lifetime — A sample of 9 x 10 6 atoms at 
800 nK confined in the horizontal dipole trap has a life- 
time of ~ 30 s. When adding the dimple beam to obtain 
settings identical to the ones used in the experiment de- 
scribed in the Letter, the lifetime is reduced to 3.4(1) s. 
We believe that this lifetime is mainly limited by 3-body 
collisions in the dimple region, as the reservoir constantly 
replenishes the dimple population. The cooling light in- 
duces an additional decay, which is strongly dependent 
on intensity and detuning. For the values used in the 
experiment, the lifetime of the reservoir is reduced to 
about 2.5 s. 

BEC lifetime — We observe that the lifetime of our 
BEC is linked to the lifetime of the reservoir gas. The 
reason is that any atom loss from the BEC is quickly and 
continuously replaced by atoms from the reservoir. A 
BEC exists as long as the total atom number in the sys- 
tem is above ~ 10 6 atoms. In order to measure the bare 
lifetime of the BEC without constant replenishment, we 
create a BEC, remove the atoms in the reservoir, and 
turn off the cooling light and the transparency light. 
The lifetime of this BEC is 1.0(1) s, probably limited 
by heating and loss resulting from inelastic 3-body colli- 
sions. The same result is obtained for an identical BEC 
created by conventional evaporation and re-compressed 
into the dimple. By contrast, the lifetime of a BEC in 
presence of the reservoir, the transparency light, and the 
cooling light corresponds to the reservoir lifetime of 2.5 s. 

Performance of the transparency beam — The 
transparency beam is absolutely necessary to protect the 
BEC from destruction by the cooling light. With the 
transparency beam turned off and the BEC subjected 
to cooling light, we measure a decay rate of about 1ms, 
and the BEC is completely destroyed in less than 3 ms. 

To demonstrate the ability of the transparency light 
to protect the BEC from resonant photons, we measure 
the lifetime of a BEC after removal of the reservoir and 
in presence or absence of both, the transparency and the 
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laser cooling light. The lifetime is the same in both cases, 
1 s. We conclude that the BEC lifetime is not limited by 
scattering of laser cooling photons if the transparency 
light is present. We give an upper bound of 1 s _1 for the 
scattering rate of laser cooling photons, which is consis- 



tent with measurements presented in Sec. IV 



The transparency beam itself has no negative effect 
on the BEC: lifetime measurements let us deduce an ab- 
sorption rate of photons from this light field to be well 
below 0.1 s _1 . 

Vacuum lifetime — Limitations to the lifetime of 
ultracold atoms in a dipole trap (beyond 3-body losses 
and residual evaporation) are off-resonant scattering of 
dipole trap photons, background gas collisions, and res- 
onant stray light. We measure the lifetime of atoms at 
low density in a deep dipole trap to be 120 s, much longer 
than any other timescale. Note that all of the lifetimes 
evaluated here are between one and three orders of mag- 
nitude larger than the BEC formation time. 



III. ABSENCE OF EVAPORATION 

The novelty of our work lies in the fact that BEC is 
achieved purely by laser cooling and thermalization be- 
tween atoms, and not by evaporation and associated loss 
to the outside environment. Standard evaporative cool- 
ing of 84 Sr has been used in past experiments and can 
be very efficient: for a factor of 10 loss in atom num- 
ber, more than three orders of magnitude in phase space 
density can be gained [57]. Thus, already little evap- 
orative atom loss can increase the phase space density 
considerably. It would require an infinitely deep trap to 
rigorously exclude evaporation. Here, we quantify the 
flux of evaporated atoms and ensure that evaporation is 
reduced to a negligible level. 

We observe an overall atom loss that leads to a lifetime 
of the sample of about 3 s. This loss is not a measure of 
evaporation, as it includes various other loss processes as 
well: 3-body recombination, loss by scattering of cooling 
photons, light-assisted collisions, and off-resonant optical 
pumping into the 3 Po,2 states by the transparency light. 
We therefore need a specific probe for evaporated atoms. 

It is also important to understand that the tempera- 
ture of atoms in the reservoir is set entirely by the cooling 
light. This is true even if efficient evaporation would be 
taking place in the reservoir: in this case, the cooling 
light would heat the sample. In the following, we will be 
concerned only about evaporation in the dimple region, 
where the sample is transparent to cooling light. 

The atoms are supported against gravity by the hori- 
zontal trap, and further confined by the dimple beam, 
which does not support the atoms against gravity. 
Atoms from the horizontal beam will preferentially evap- 
orate downwards, aided by gravity. If they do so within 
the dimple region, they are attracted by and guided in 
the dimple beam. This beam provides superb confine- 
ment of evaporated atoms and enables us to detect even 
minute atom numbers in absorption imaging. The dim- 



ple beam has a potential depth of ks x 2.6 fj,K and there- 
fore selectively guides only low-energy atoms from evap- 
oration, but not atoms originating from inelastic colli- 
sions. 

The amount of atoms leaking from the dimple region 
into the dimple beam will depend on the sample temper- 
ature, the power of the horizontal beam, and the power 
of the dimple beam. We set the temperature to the one 
used in the experiment, vary the power of the horizon- 
tal and the dimple beam independently, and measure 
the atom leakage. The result of such a measurement is 
shown in Fig. [S4j The data points denote combinations 
of parameters where the atom flux is about 50 atoms/ms, 
which is our detection threshold. The straight line is 
a linear fit to the data. Note that for a leakage of 
50 atoms/ms and a BEC formation time below 60 ms, 
less than 3 000 out of 9 x 10 6 atoms evaporate in the dim- 
ple region. For the experiments described in the Letter 
we verify that no discernable evaporation takes place at 
any time during the optical cooling. 

We explore the entire parameter region shown in 



Fig. S4 and find no dramatic difference in BEC formation 
time or final BEC atom number, again indicating that 
any residual evaporation has no measurable contribution 
to the BEC formation. 



IV. SHIFTING ATOMIC AND MOLECULAR 
TRANSITIONS 



A. Introduction 

We intend to cool atoms in the dipole trap with light 
resonant to the 1 So — 3 Pi transition, but only in the 
region outside the dimple. To do so, we render the atoms 
in the dimple region transparent for the cooling light by 
locally shifting the 3 Pi state out of resonance. Cooling 
of atoms in the reservoir is performed on the red side of 
the 1 Sq — 3 Pi transition, so we need to locally shift the 
3 Pi state towards higher energies by at least a few 10 
linewidths. This is done with the transparency beam, 
which spatially overlaps with the dimple region. 

To characterize the system, we work in similar condi- 
tions as used in the experiment described in the Letter. 
Atoms are tightly confined in a crossed-beam dipole trap, 
consisting of the horizontal trap and the dimple, and we 
evaporate into a pure BEC to reach both high densities 
(to enhance possible photoassociative loss) and low trap 
depths (such that atoms can leave the trap with very 
few photon recoils). We then illuminate the atoms with 
circularly polarized cooling light, using only an upward 
propagating beam. The intensity of the cooling light is 
1.2 / sat and the duration of illumination is 10 ms. A mag- 
netic field of 90 mG magnitude and vertical orientation 
is present. A typical loss spectrum is shown in Fig 



Sr, 



note the broad photoassociation (PA) feature, which cor- 
responds to the v = — 1 state of the 1(0+) potential [51] . 
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Figure S4: Detection of residual evaporation. Hot atoms from the dimple region might leave the dipole trap along gravity. 
When doing so, they are funneled into the dimple beam and can be detected by absorption imaging (insets). The amount 
of evaporation depends on the power in the horizontal and dimple trap beams, illustrated by the three insets showing atom 
fluxes of 350, 120, and < 50atoms/ms. Our detection limit is about 50atoms/ms, and conditions with atom fluxes that touch 
this limit are indicated by the black data points. The straight line is a linear fit to the data. Experiments are performed in 
the region below this line, indicated by the star. The images are averages of 10 experimental realizations, and the color scale 
is adjusted to maximize contrast of small signals. 




sition. The light is linearly polarized, propagating at an 
angle of 15° to vertical. The waist is set to 35 (im to en- 
sure a uniform light intensity across the sample. We take 
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Figure S5: Loss spectrum of atoms confined in a dipole 
trap and illuminated by a + -polarized cooling light. A small 
guiding field of 90 mG splits the Pi state into three clearly 
visible m j states. The broad feature around —400 kHz is due 
to photoassociation (PA). The visibility of ir and a~ transi- 
tions indicates that the cooling light is not perfectly circularly 
polarized. In the experiment described in the Letter, cooling 
is performed on the red side of the a + transition, indicated 
by the arrow. The line is a guide to the eye. Note that the 
detuning here is calculated from the resonance position at 
zero B-field, which is determined to better than 10 kHz. 



B. Shifting of resonances 

In a next step, we illuminate the atoms with trans- 
parency light tuned 15 GHz blue to the 3 Pi — 3 Si tran- 



loss spectra as the one shown in Fig.[S5]ior various inten- 
sities and observe that the four loss features experience 
different Stark shifts and move differently with intensity; 
see Fig. |S6| The 7r-transition moves rapidly towards blue 
detuning, whereas the a + and PA lines approach asymp- 
totic values less than 100 kHz from resonance. 

For our cooling strategy, we intend to work at a detun- 
ing of about —150 kHz from the resonance frequency at 
zero B-field. It seems possible to move the a + transition 
a few 100 kHz away from this frequency. Unfortunately, 
the broad PA transition moves into exactly this region, 
potentially causing undesired losses given the high den- 
sity in the dimple region. The transitions broaden con- 
siderably with increased transparency beam intensity. 



C. Polarization dependence 

So far, we have not considered the polarization of the 
transparency beam, which was linear in the previous 
measurement. This beam is almost co-aligned with the 
magnetic quantization axis, such that its field can be de- 
composed into equal amounts of <r + and er~ -light (but no 
7r-light) in the reference frame of the atoms. The three 
mj states of the 3 Pi manifold, as well as the PA line 
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Figure S6: Resonance positions in dependence of transparency light intensity. Spectra as the one shown in Fig. |S5| are taken 
for various intensities of the transparency beam, and the resonance positions recorded. The a~ transition is too weak to 
be tracked. Inset: a magnification of the region around the cool ing frequency, denoted by the dotted line. An intensity of 
0.6 MW/m 2 (arrow) is used to take the data presented in Fig. S7 



tied to the m 



,/ = +1 state, couple differently to the 3 Si 
states; see Fig. S7'a). This helps to understand the dif- 
ferent shifts of the it- and cr + -transitions in Fig. S6 and 
we expect the shifting of the individual states to strongly 
depend on the polarization of the transparency light. 

In a next step, we fix the intensity of the transparency 
light intensity to 0.6 MW/m 2 and take scans at vari- 
ous polarizations. The results are shown in Figs. |S7[ b) 
and (c). The loss features change dramatically in po- 
sition, width, and amplitude depending on the polar- 
ization of the transparency beam. The er + -transition 
and its associated PA line are maximally shifted by a 
cr~-polarized transparency beam. This can easily be un- 
derstood from the fact that c + -polarized light cannot 
couple the mj = +1 state to the 3 Si manifold. We find 
that already small variations in the polarization decrease 
the performance considerably. Note that in the experi- 
ment described in the Letter, we use an intensity of the 
transparency light that is larger by a factor of about 10. 

We will now focus on the region of interest around 
the cooling frequency and highlight the importance of 
the polarization of the transparency light under condi- 
tions close to the ones used in the experiment described 
in the Letter. We increase the intensity of the trans- 
parency beam to 7 MW/m 2 and make the cooling pulse 
25-times stronger (5-times longer and 5-times more in- 
tense) to amplify even small losses. Two spectra for 
cr~ and near-linear polarization are shown in Fig. S7 d). 
Not only does the loss feature move to positive detun- 
ing with increasing circular polarization, also its width 



and amplitude decrease dramatically. The small remain- 
ing loss feature is about 250 kHz = 35 T away from the 
desired cooling frequency. We speculate that this is the 
mj = — 1 state, which is unaffected by transparency light 
of o~ -polarization and driven by small er~-polarization 
components in the cooling light. 



D. Conclusion 



We have shown that the magnetic substates of the 
3 Pi state, as well as the highest molecular level of the 
1(0+) potential, can be shifted by a light field slightly 
detuned from the 3 -Pi — 3 «5i transition. In particular, 
we find that the three Zeeman lines and the first PA 
line move very differently in dependence of the polariza- 
tion of the transparency laser. We find that all lines can 
be moved at least 250 kHz away from the desired cool- 
ing transition, reducing the photon scattering rate by at 
least three orders of magnitude. The strongest lines are 
even moved by a few MHz. This transparency beam al- 
lows us to cool atoms in the reservoir, while keeping the 
photon scattering rate of atoms in the dimple below the 
rate of 3-body collisions, which eventually determine the 
lifetime of the BEC. Further investigations are needed to 
gain a quantitative understanding of the AC Stark shifts 
in dependence on transparency laser frequency, intensity, 
and polarization. 
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Figure S7: Resonance positions in dependence of the polarization of the transparency light, (a) The three mj states of 
the 3 Pi manifold couple differently to light on the 3 Pi — 3 £i transition (green arrows). The mj — +1 state, for example, 
does not couple at all to <r + -polarized light due to the absence of any mj = +2 state in the 3 Si manifold. The cooling light 
(red arrow) is always <r + -polarized, with small admixtures of other components due to experimental imperfections, (b) A 
loss spectrum of atoms illuminated by the transparency light and subjected to a pulse of cooling light. The frequency of the 
cooling light is scanned, and the number of remaining atoms is recorded. The three loss features correspond to the PA-, o + -, 
and 7r-transitions. This spectrum is taken with a near-linear polarization of the transparency light and corresponds to the 



data set indicated by an arrow in Fig. S6 (c) We take spectra at various polarizations, covering the complete range between 
<r + -circular (bottom) to linear and cr _ -circular polarization (top). The fraction of remaining atoms is encoded in the color, 
with red color denoting high atom loss. The dotted white line indicates the frequency of the cooling light, (d) Loss spectra for 
two extreme values of the polarization, taken with the intensity of the transparency light used in the experiment described in 
the Letter, and applying a 25-times stronger cooling pulse compared to the measurements shown in (b) and (c). The spectra 
taken with intermediate polarizations (not shown) smoothly interconnect between these two curves. The cooling frequency is 
again indicated by an arrow. 



V. DENSITY DISTRIBUTION MODELS 

We now determine the density distributions of our 
samples by a numerical model and by a simple ana- 
lytical model. Both models assume that the gas is in 
thermal equilibrium and neglect the influence of laser 
cooling light on the density distribution in the reservoir. 
The models are compared to the experiment and used 
to determine important quantities, such as the elastic 
collision rate. 



A. Self-consistent mean-field equations 

We use a self-consistent mean-field model to describe 
the density distribution of the BEC n$(v) and of the 
thermal gas n th (r). The thermal gas is described as a 
gas of non-interacting particles residing in the potential 
^th(r) = U ex t(r)+2g(n (r)+n t h(r)), which is the sum of 
the external trapping potential U cx t and the mean-field 



potential of all atoms. Here g = 47rfi 2 a/m, where a is 
the scattering length and m the mass of the atoms. The 
potential C/ GX t is assumed to be zero at its deepest point. 
In the semiclassical approximation, the thermal density 
distribution in dependence of the potential is given by 

where Ath = hj s/^hvwM^T is the thermal de Broglie 
wavelength and fi the chemical potential. The func- 
tion (z) is a polylogarithm function, g a {z) = 
S^Li z n /n a - The spatial density distribution is then 

nth(r) = n th ([/t h (r)). (SI) 

The BEC distribution is obtained in the Thomas-Fermi 
(TF) approximation 

7i (r)=max[( M -LVr))/s,0], (S2) 
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where the potential Uo(r) — U ex t(r) + 2gn t h(r) is the 
sum of external potential and the mean-field potential 
of the thermal atoms. 



B. Numerical model 



To determine the density distributions no and nth, we 



self-consistently solve Eqns. SI and S2 The solution is 
obtained by a numerical calculation for given values of 
the chemical potential /x and the temperature T . The 
calculation starts from initial density distributions that 
are zero everywhere, and obtains better and better ap- 
proximations for the distributions by iterating Eqns. |S1| 
and |S2| The potential is described by a sum of the 
gravitational potential and the potentials created by the 
dipole trapping beams, which we assume to be Gaussian. 

To speed up calculations, we first calculate tables con- 
taining no(U) and n t h(U) for a range of relevant poten- 
tial depths U. We then interpolate these tables to obtain 
no(U) and n t h(U) for any U. The spatial density distri- 
butions no(r) and n t h(r) are obtained in the local den- 
sity approximation using realistic potentials U(r). To 
avoid numerical artifacts, the self-consistent solution of 
Eqns. [ST] and [S2] is obtained in two stages. During the 
first stage, the mean-field of the thermal atoms is not 
included in the potential that the thermal atoms experi- 
ence. After 20 rounds of iteratively solving Eqns. [Si] and 
|S2[ we linearly ramp on this potential contribution in a 
second stage consisting of 25 rounds. In every round, the 
minimum f7 t ™ m of [/th is determined. Since the presence 
of a BEC is assumed, n t h(U™ n ) equals n c which is used 
as a boundary condition in the determination of n t h(i"). 
To avoid oscillations in the calculation, we low-pass fil- 
ter U^ ln and use the filtered value in the calculation of 
rcth(r). 

In the experiment, we determine the atom number in 
the reservoir N ICS and the atom number in the dim- 
ple iVdimpie by performing a double-Gauss fit to the 
density distribution recorded after a short expansion 
time. To determine similar values from the numerical 
calculation, we define the reservoir density distribution 

^th,rcs 

(r) = n t h (Kes (r)), where U ICS (r) is a potential 
without the contribution of the dimple, but the same 
potential values outside the dimple region. The atom 
numbers in BEC and reservoir are determined by in- 
tegrating rio(r) and nth,res( r ), respectively. The ther- 
mal atom number in the dimple is obtained by inte- 
grating n t h,dim P ic(r) = nth(r) - n t h,res(r). The radii Ri 
(i E {x, y, z}) of the BEC are defined by the locations at 
which the BEC density drops to zero. 

An important quantity which we want to extract from 
our model is the elastic collision rate of the the rmal gas. 
It is given by r e i itri (^) = nvcr, where v = A^ksT /nm 
and a = 87ra 2 . The peak collision rate of the thermal gas 
Tel, th, peak is reached at the edge of the BEC, where the 
gas has the critical density. The timescale of thermal- 
ization between the gas in the dimple and the reservoir 
gas is limited by the elastic collision rate in the reservoir, 



outside the dimple region. This rate is highest just to the 
side of the dimple in the xy-plane, and determined from 
n t h at that location, which we denote as n t h,edgc- Assum- 
ing a Gaussian reservoir density profile in the ^-direction 
and a homogeneous distribution in the xy-plane, the av- 
erage elastic collision rate in the reservoir close to the 
dimple is r el;th>res = n th ^ cdgc va/\/2. 



C. Analytical model 

Before we present the results of the numerical model, 
we discuss a simple analytical model. We neglect the 
mean-field of the thermal atoms and approximate the 
dimple potential by a harmonic potential with trap os- 
cillation frequencies /j. The BEC resides entirely in the 
dimple and has the usual TF inverted parabola shape. 
The chemical potential is determined from the BEC 
atom number No by 



Tko ( l5N a 

Oho 



2/5 



where = ^fi/mw and u = 27r(f x f y f z ) 1 / 3 . The 
radii of the BEC are given by the TF radii i?TF,j = 
\Jl\ilm{2-n]i)' 1 and the peak density is no(0) = ji/g. 

Because of different correlations in the BEC and the 
thermal gas, the BEC mean-field is twice as strong for 
thermal atoms than for atoms in the BEC. Thermal 
atoms in the center of the trap experience a "Mexican 
hat" potential, created by the external confinement com- 
bined with the BEC mean-field. In the center of the 
BEC, this mean-field reaches 2</rio(0) = 2/i. For our typ- 
ical parameters (see Tab. SIII), the chemical potential is 
on the order of half the temperature, /i ~ 0.5fcsT, which 
means that the thermal distribution is strongly modified 
by the BEC mean-field. The density of the thermal gas 
reaches the critical density n c — n t h(/i) = 2.612/Aj^ at 
the edge of the BEC. The dimple potential has a finite 
size in the xy-plane and a finite depth U*j?^ lc . At the 
edge of the dimple, the density is n t h,cd g c = n t h(U^ plc ). 
The elastic collision rate in the reservoir is given by 
Tel, th, res = r cl ^h(«th,odgc) and the peak elastic collision 
rate of the thermal gas is again r o i !t h, P cak = r e i )t h(n c ). 

We can estimate the atom number in the reservoir 
Arcs by describing the reservoir as a harmonic poten- 
tial with trap oscillation frequencies / ros ,i- The cen- 
tral density of the reservoir is about n t h,cdgc, which is 
much lower than the critical density. Therefore the gas 
in the reservoir is well described by a classical ther- 
mal distribution. With this approximation, we ob- 
tain N ICS = nth,cdge(27rfc B r/m) 3 / 2 /a;j 5 0S , where w ros = 

2'K'(frcs,xfrcs,yfvcs,z) ^ • 



D. Comparison to the experiment 

We now compare the experiment and the numerical 
model. Since the model can only describe a thermalized 
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Figure S8: Comparison of experiment (blue data points) and numerical model (solid lines). Shown are BEC atom number 
(a) and temperature (b) in dependence of the total atom number. Experiment and calculation agree well. A peculiar feature 
is that the BEC atom number is nearly proportional to the total atom number. An important contribution to this behavior 
is the increase in temperature with the atom number observed in the experiment. For comparison, a calculation for constant 
temperature is shown as dashed line. Here, the BEC atom number increases much more steeply with total atom number. 
Results of the models matched to the star-shaped data point are given in Tab. |SIlTl and Fig. |S9| 
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Figure S9: Potential (a, b) and density (c, d) along the y-direction determined by the numerical model, using parameters 



that correspond to the star-shaped data point of Fig. S8 In (a) and (b), the reservoir potential is shown as a blue, dashed 
line and the full potential, including the dimple, as a red line. Panels (c) and (d) show the density distributions of the thermal 
gas n t h (red line) and of the BEC no (black line). The dashed blue line shows the density of thermal atoms in absence of the 
dimple potential n t h,ro S . 
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sample without the influence of the laser cooling beams, 
we perform experiments corresponding to this situation. 
The experimental sequence is identical to the one de- 
scribed in the Letter, up to the moment before the dim- 
ple is ramped on. Here, we switch off laser cooling and 
we do not use the transparency laser. After ramping on 
the dimple, we wait 200 ms to let the system thermalize. 
Then we record as usual time-of-flight absorption im- 
ages with 24 ms expansion time, from which we extract 
temperature, total atom number and BEC atom number. 
We perform this measurement for different loading times 
of the metastable state reservoir at the beginning of the 
experimental sequence. A longer loading time leads to a 
higher atom number and, as a consequence of the laser 
cooling process, to a higher temperature. The result of 
the experiment is shown in Fig. |S8| In the examined 
range from 3.5 x 10 6 to 8 x 10 6 atoms, the temperature 
changes from 880 nK to 990 nK. The atom number in the 
condensate is nearly proportional to the total number of 
atoms iVtot and is well described by N = 0.027V t ot- 

The numerical model needs the potential as input. To 
describe the potential, we sum the contributions of grav- 
ity and of the potentials of the three Gaussian optical 
dipole trap beams used in this experiment. The pa- 
rameters of the beams, determined from trap oscillation 
and beam power measurements, are given in Tab. |SII| 
The precision with which we know the potential shape 
is not quite high enough to give a good match between 
the model and the experiment. We therefore vary one 
potential parameter, here the dimple depth, and match 
the model to the experiment around iVtot = 6 x 10 6 . 
Experiment and model agree for a dimple depth that 
is 10% larger than the measured one. A good match 
between model and experiment could have alternatively 
been achieved by varying other parameters or by cor- 
recting a possible systematic error in the experimental 
data. A candidate systematic measurement error, which 
would be sufficient to bring the numerical model and the 
experiment into agreement, is an overestimation of the 
temperature by 5%. This hypothetical systematic shift 
is smaller than our experimental uncertainty. 

After matching the model to the data in a single point, 
we can now explore the whole range of data. The model 
requires as input parameters the chemical potential and 
the temperature. For two data points near 3.5 x 10 6 
and 8 x 10 6 atoms, we determine the chemical potential 
that corresponds to the measured BEC atom number. 
This procedure results automatically in a good match 
for the total atom number, which already shows that the 



Table SII: Parameters of the dipole trap beams used for the 
experiments described in this section. 



beam 


waist x 

(H 


waist y 

w 


waist 2 

w 


P 

(mW) 


U/k B 
( M K) 


reservoir 




298.2 


17.32 


1748 


11.73 


vertical 


297 


297 




743 


0.29 


dimple 


22.41 


22.41 




39.6 


2.73 



model describes the data well. Instead of determining 
the chemical potential for each measurement point, we 
now simply linearly interpolate temperature and chemi- 
cal potential between these two extreme points and cal- 
culate No and iV to t for many values in between. The re- 



sult is the solid line in Fig. S8 and shows a good overall 
match between experiment and calculation. In a second 
calculation, we keep the temperature always at 930 nK 
and vary only the chemical potential; see the dashed line 
in Fig. |S8| The mismatch between calculation and data 
shows the importance of the small temperature change 
observed in the experiment between low and high atom 
number. 

We now examine the data point marked by a star in 
Fig. [S8] in more detail. We chose this data point be- 
cause the atom number in the BEC is comparable to the 
one achieved in the experiments discussed in the Let- 
ter. As input for the models, we use the measured tem- 
perature and BEC atom number and adapt the chemi- 



cal potential accordingly. Figure [S9| shows the potential 
and the density profiles of BEC and thermal gas along 
the x-direction as determined by the numerical calcu- 
lation. The BEC mean-field clearly has a strong effect 
on the thermal distribution, pushing thermal atoms out- 
wards. Table [SlTT| compares experimental data with the 
two models and shows a general good agreement. 

Some information extracted from the models is worth 
further discussion. We expect that the conditions in the 
experiments presented in the Letter are comparable to 
the experiments presented here within the volume illumi- 
nated by the transparency laser. Only outside this vol- 
ume, laser cooling will change the conditions and make 
a comparison more complicated. An indication for this 
is that the atom number in the dimple, which was not 
measured here, is comparable to the one observed in the 
experiment of the Letter. 



Table SIII: Comparison of the experiment and the models. 
The values given correspond to the star-shaped data point 
in Fig. |S8| The trap oscillation frequencies of dimple and 
reservoir are f x = 244 Hz, f y = 247 Hz, f z = 624 Hz, f tes>x = 
8.4 Hz /res,!/ ~ 37 Hz, and / re s,z = 622 Hz. 



quantity 


experiment 


numerical 


analytical 


T 


944(50) nK 


944 nK 


944 nK 


n/k B 




478 nK 


405 nK 


-/Vbec 


111 x 10 3 


111 x 10 3 


111 x 10 3 


dimple 




654 x 10 3 


2.6 x 10 B 


N reservoir 




5.0 x 10 6 


•Wtotal 


5.8 x 10 6 


5.8 x 10 6 




no(0) 




1.0 x 10 15 cm -3 


8.6 x 10 14 cm 


n c 




3.5 x 10 14 cm~ 3 


3.5 x 10 14 cm 


^th.cdgc 




9.9 x 10 12 cm -3 


8.7 x 10 12 cm 


V 




21.8 /im/ms 


21.8 (jm/ms 
2700 s" 1 


lei, th, peak 




2700 s- 1 


Tel, th, res 




50 s" 1 


50 s" 1 


RtF.x 




5.1 fim 


5.8 fim 


Rtf, v 




5.0 /im 


5.8 fim 


RtF.z 




1.9 /im 


2.3 /im 



-3 
-3 
-3 
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Since the reservoir contains nearly an order of mag- 
nitude more atoms than the dimple, the density distri- 
bution of the gas in the reservoir and the temperature 
of the system will barely change when ramping the dim- 
ple on. Therefore the ratio of n. t h,cdgc/^c, which is 35, 
gives an idea of the thermal gas density increase and the 
phase-space density increase provided by the dimple. 



Finally, we note that the average elastic collision rate 
in the reservoir close to the dimple corresponds to about 
five collisions during the ~ 100 ms required for thermal- 
ization. This number compares well with the ~ 3 col- 
lisions required to obtain thermal equilibrium as deter- 
mined from Monte Carlo simulations. 
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